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Abstract 



The problem of front propagation in flowing media is addressed for laminar 
velocity fields in two dimensions. Three representative cases are discussed: 
stationary cellular flow, stationary shear flow, and percolating flow. Produc- 
tion terms of Fisher-Kolmogorov-Petrovskii-Piskunov type and of Arrhenius 
type are considered under the assumption of no feedback of the concentration 
on the velocity. Numerical simulations of advection-reaction-diffusion equa- 
tions have been performed by an algorithm based on discrete-time maps. The 
results show a generic enhancement of the speed of front propagation by the 
underlying flow. For small molecular diffusivity, the front speed Vf depends 
on the typical flow velocity U as a power law with an exponent depending on 
the topological properties of the flow, and on the ratio of reactive and advec- 
tive time-scales. For open-streamline flows we find always Vf ~ U, whereas 
for cellular flows we observe Vf ~ C/ 1 / 4 for fast advection, and Vf ~ C/ 3 / 4 for 
slow advection. 
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I. INTRODUCTION 



Interface motion and front propagation in fluids occur in many different areas of interest to 
science and technology. Among the most important examples we mention chemical reac- 
tion fronts in liquids, population dynamics of ecological communities (e.g. plankton in the 
ocean), atmospheric chemistry (ozone hole) and flame propagation in gases M. 
The mathematical description of those phenomena is based on partial differential equations 
(PDE) for the evolution of the concentration of the reacting species and the evolution of 
the velocity field ||. In principle the two PDEs (for the reactants and the velocity field) 
are coupled, often in a nontrivial way. An example is given by a recent study of reactants 
coupled to the Navier-Stokes equation by a Boussinesq term ||. A mathematical simplifi- 
cation can be obtained assuming that the reactants do not influence the velocity field which 
evolves independently. In such a limit the dynamics is still completely nontrivial and it is 
described by a so-called advection-reaction-diffusion equation. In the most compact model 
one considers a single scalar field #(x, t), which represents the fractional concentration of 
products. The field 9 has a zero value in the regions containing fresh material only, 9 is unity 
where the reaction is over and there are only inert products left. In the region where the 
production takes place and reactants and products co-exist, the field 9 assumes intermediate 
values. 

The evolution of 9 in a reacting fluid with molecular diffusivity Dq is described by the PDE 

l e + (u-V)9 = D V 2 9 + ±f(9) (1) 

where u(x, t) is an incompressible velocity field. The second term on the r.h.s. of Eq. ([[]) 
describes the production process, characterized by a typical time r r . The shape of f(9) 
depends on the phenomenon under investigation, and we will consider two relevant functional 
forms: 

1) f(9) = 9(1 — 9), or more generally any convex function (f"(9) < 0) such that 
/(0) = /(l) = 0, /'(0) > and f'(l) < 0. This is called the Fisher-Kolmogorov- 
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Petrovskii-Piskunov (FKPP) nonlinearity |4]]. The production term is proportional to 
the concentration of reactants, 1 — 9, and to the concentration of products, 9. 

2) f{9) = e~ dc l d {l — 9). This is the Arrhenius nonlinearity ||. In this case, the solution 
9 = is only marginally unstable, and the parameter 9 C plays the role of an activation 
concentration, since practically no production takes place for products concentrations 
below that threshold. Production is still proportional to the concentration of fresh 
material. 

We will always consider initial conditions such that #(x, 0) — > 1 exponentially fast as x, the 
horizontal coordinate, approaches — oo, and #(x, 0) — > exponentially fast as x — > oo. The 
initial profile of 9 has no variation along the transversal axis. The choice of such an initial 
condition, for which the concentration of products has a noncompact support, is of interest 
because it suppresses all possible flame-quenching effects that may appear in case 2) (see 
e.g. ||). This allows a direct comparison between FKPP and Arrhenius production terms. 
There exists a huge literature for the case u = (see e.g. ||). In that case the physical 
mechanism for front propagation resides in the the combined effect of diffusion and produc- 
tion. Let us indeed consider a one-dimensional situation where a reservoir of fresh material 
is located on the right side, whereas on the opposite side we have only inert products. At 
the boundary between the two phases, diffusion mixes fresh material and inert products, 
broadening the interface. Then, production raises the level of the concentration of products, 
thus shifting the interface, to the right in this case. The final result is a front propagating 
from left to right, eating out fresh material to leave behind digested, inert products. The 
front speed at large times reaches a limiting value Vq. For the FKPP nonlinearity one has 
the exact result 




(2) 



whereas for a generic f(9) one has the bounds 
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It has to be remarked that the convergence to the limiting velocity is extremely slow for 
FKPP production |||9[|, therefore this case requires special attention especially for nonuni- 
form flow. 

In the presence of a moving medium, i.e. u / 0, one expects that the front propagates 
with an average limiting speed Vf. A problem of primary interest is to determine the de- 
pendence of Vf on the properties of the velocity field u [ To|,ll| . In this article we consider 



front propagation in simple laminar flows (shear flow and systems with cellular structures) 
which, in spite of their apparent simplicity, show intriguing behaviour [12|-|l5l. For a given 



structure of the flow field, we aim to explore the dependence of Vf on relevant parameters, 
such as the typical flow velocity U and the production time-scale r r . In terms of adimen- 
sional quantities, we look for an expression for the speed enhancement Vf/Vo in terms of 
the Damkohler number, Da = L/(Ur r ), which measures the ratio of advective to reactive 
time-scales, and in terms of the Peclet number, Pe = UL/Dq which expresses the relative 
weight of advection and diffusion. We will mainly be interested in the case of large Pe, to 
highlight the combined effects of advection and reaction. The two regimes Da <C 1 and 
Da ^> 1 are quite different in nature: in the first case, typical for slow reaction rates or fast 
advection, the front interface is distributed over several length-scales L, and for this reason 
it goes under the name of "distributed reaction zone" regime ||16|| ; in the second case, the 



front is thin compared to L, and it propagates according to a Huyghens-like principle, hence 



the name of "geometrical optics" regime [16]. We will provide a detailed analysis of these 



two regimes, highlighting their differences. In the context of the thin front regime we can 
mention, among the many contributions, the works on the G-equation approximation and 



its relation with the "geometrical optics" regime [19|]2(| , the work on turbulent flows ||10| . |11 
and the numerical study of front propagation in synthetic turbulence [f23[|. 
A hint to the effect of an underlying flow on front propagation is given by the observation 
that the front speed increases with the square-root of molecular diffusivity. It is well-known 
that diffusive transport is always enhanced by incompressible flow, resulting in an effective 



diffusion coefficient _D e fr > -Do ||27| , |32| , |33| ■ Therefore, it is reasonable to expect that the front 



speed will be enhanced, too. This physical argument can be upgraded to a mathematically 
rigorous statement in the case of a "slow" reaction, that is for Da -C 1 (see Section [TT]). 
From the mathematical viewpoint, there exist lower bounds to the speed of the front which 
confirm the expectation that the flow enhances front propagation p~7|JT8|] . These bounds 
take different forms according to the topological structure of the flow field. One can dis- 
tinguish two main classes: cellular flows, characterized by having closed streamlines only, 
and percolating flows, which possess open streamlines (shear flows are a particular case of 
this second class, having only open streamlines). For cellular flows, it has been shown that 
Vf/Vo > C x Dar x f 2 + C 2 for Da > 1 and V f /V > C x Da- x l h + C 2 for Da < 1, with constants 
Ci,C 2 depending on the shape of the production term. For percolating flows, the lower 
bound is expressed as Vf > K-JJ where again K x is a constant depending on /. The most 
general upper bound, valid for flows of both classes, is Vf < Vq + K 2 U |17|]. We therefore 



see that percolating flows are constrained to a linear dependence on the stirring intensity 
U (see Section [IV A| for numerical results). Cellular flows have more intriguing properties, 



as we will see in Section |1V BE for fast advection, Da <C 1, the front speed depends on the 
flow velocity as Vf oc f/ 1 / 4 , to be compared with the lower bound prediction oc U 1 ^ and the 
upper bound oc U; for slow advection, Da ^> 1, we obtain Vf oc U 3 ^ to compare with the 
lower bound oc U 1 ^ 2 and the upper bound oc U. It is clear that these bounds do not provide 
a sharp evaluation of the front speed in the case of a cellular flow. 

To close the overview of mathematical results, we anticipate that a different upper bound 
for the front speed enhancement can be obtained by reformulating the solutions of Eq. (|]) 
in terms of a path-integral (see Section [II] for details, and references therein). This bound 
yields an expression similar to the one obtained in the absence of any flow, but with an 
effective diffusion coefficient D e g replacing the molecular one. Explicitly, we show that 



T r . „ D eS f(6) 

V f <2^ — sup— (4) 

where the dependence of _D e fr on the flow parameters and on the molecular diffusion can 
be derived by the analysis of Eq. ([J) when the production term has been switched off. 



For a cellular flow, we have the result D e s ~ y/ULDo, fl4]Jl3f| , whereas for a shear flow 
D cS ~ (UL) 2 /Dq 0. Inserting these latter expressions in Eq. we obtain the behaviour 
V} oc t/ 1 / 4 for the cellular flow also obtained by |36[] and Vf ocU for the shear flow, re- 



markably close to the observed ones for Da -C 1. Furthermore, the upper bound @ is sharp 
in the regime where homogenization techniques apply (see Section [TI| and Appendix |C|). In 
other words, for fast advection the effect of the underlying flow can be compactly expressed 
in the renormalization of the diffusion coefficient. On the contrary, for slow advection in a 
cellular flow, the front speed departs significantly from the upper bound, with an increase 



in front speed less prominent than in the fast advection regime. We will see in Section [IV B 
that the main physical mechanism accounting for this depletion is the appearance of an 
effective reaction term as a consequence of the joint effect of advection and reaction. 
These observations lead us to argue that the effect of a stirring velocity on front propagation 
can be in general summarized in the renormalization of the two relevant parameters: (i) an 
effective diffusivity, which is always larger than the molecular one for an incompressible 
flow, and (ii) an effective production term, which is slower than the microscopic one for slow 
advection. 

The remainder of the paper is organized as follows: in Section [TI|, we derive some general 
results valid for all advecting flows; Section |lj presents the algorithm employed for the 
numerical solution of Equation ([I]); in Section |TV] we discuss the results of numerical simu- 
lations; Section [V| is devoted to conclusions and discussion. Technical and numerical details 
are treated in the Appendices. 

II. UPPER BOUNDS TO FRONT SPEED 

In this Section we show how to establish the upper bound (|]) for the speed of front propaga- 
tion in a generic incompressible flow and a generic production term. This result is the conse- 
quence of the deeply rooted link existing between front propagation and advective transport. 
In other words, we will exploit the relationship of the solutions of Eq. ([[]) to the solutions of 
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the same equations in absence of production terms. This will yield the constraint (^) involv- 
ing the front speed Vf, the effective diffusion coefficient D e g and the production time-scale 
r r . In general this relation is an inequality, and not a sharp functional relation. In a graphi- 
cal form, this amounts to say that "front propagation^ (advection+diffusion)+production". 
A significant exception to this general rule is given by the limit of very slow reaction (or 
very fast advection), where the bound (Q) becomes sharp. In this case, homogenization 
techniques, detailed in Appendix [A], allow to show that the problem of the determination 
of the front speed reduces to the problem of determining D c r. This is essentially due to the 
large separation of typical time-scales. 

We start our proof of the inequality (^) by recalling the fundamental relation among the 
solution of the PDE ([!]) and the trajectories of particles advected by a velocity field u(x, t) 
and subject to molecular diffusion [p5[ , f26H 

0(x, t) = (0(r(O), 0) e-tf °(*W»)."))<^ (5) 

where 

is the growth rate of 9. The average is performed over the trajectories evolving according 
to the Langevin equation 

dr(t) 
dt 



v(T(t),t) + y/2D rj(t) (6) 



with final conditions r(t) = x. The white noise term y/2Dorj(t) accounts for molecular 
diffusion. 

Since the growth rate is always bounded from above, c{9) < c max = sup^ c(9), Eq. (0) yields 
the inequality 

0(t,x) < (fl(0,r(0))exp(c max t) . (7) 

For FKPP production terms the maximum occurs for 9 = 0, that is we have c{9) < 
c(0) = f'{0)/r r = l/r r , therefore c max = l/r r . In the inequality (0), the term in an- 



gular brackets denotes the probability that the trajectory ending at x were initially lo- 
cated at the left of the front interface. Under very broad conditions, i.e. nonzero molec- 
ular diffusivity and finite variance of the velocity vector potential p?tp2| , p3[ , it is pos- 



sible to show that asymptotically the particles undergo a normal diffusion process with 
an effective diffusion coefficient _D e ff, always larger than the molecular value D . This 



is the rigorous and most general version of a statement originally due to Taylor p4 
The issue of single particle diffusion, and the problem of finding the effective diffusiv- 
ity, given a velocity field and a molecular diffusivity, has been the subject matter of a 
huge amount of work (see e.g. [^] for a recent review). In the presence of an asymp- 
totic normal diffusion, we can substitute the term (#(0,r(0)), with the gaussian result 
1 — |erfc(— xj y/2D cS t) ~ exp[— x 2 / (4D cS t)]/ ^/2nlJ^t, where the latter approximation holds 
with exponential accuracy. We thus obtain 6*(t,x) < exp [c max t — x 2 / (4_D eff t)] / \/2nD c{i t. 
It is therefore clear that at the point x the field 9 is exponentially small until a time t 
of the order of x/\/4:D cS c max . We therefore obtain the upper bound for the front velocity 



Vf < v / 4-D c frc m ax, as anticipated in Eq. ([|). 

The analytic determination of the effective diffusivity from the knowledge of the advecting 
field and the value of the molecular diffusivity is - in general - a daunting task. Nevertheless 
there is an exact result valid for all flows in the form of an upper bound for the effective 
diffusivity D e g < D (l + aPe 2 ), where a is a numerical constant that depends on the details 
of the flow p7| , |32| , |33| . Plugging this relation into @ we can derive a general upper bound 



Vf/Vo < Vl + aPe 2 (8) 
where any dependence on the flow details has been summarized in the numerical constant 



a. For large Pe we recover the bound Vf < const.?/ already discussed in Ref. |T7|]. In the 
limit of small Pe, i.e. for small stirring intensity U, this bound is in agreement with the 
Clavin- Williams relation (Vf — Vq)/Vq = (U/Vq) 2 ||10|| . On the contrary, in the same limit, 
an asymptotic behavior like (Vf — V )/V ~ (£//Vo) 4//3 , proposed in Ref. |TT|] , is ruled out 
since the rate of convergence to Vo for vanishing U has to be faster or equal to U 2 in order 
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to fulfill the bound <g). 

As anticipated above, there is a situation where the bound (|j) becomes sharp, and that is 
the limit of very slow reaction. It is easy to understand the physical reasons for this effect. 
If r r is the slowest time-scale under consideration, advection and molecular diffusion act 
jointly to build an effective diffusion process, unaffected by reaction. Diffusion decreases 
the value of concentration to lower levels before the onset of production, which then takes 
place at the maximal rate (e.g., for FKPP, at 6 ~ 0). In the limit of very slow reaction, 
basically one has that Eq. ([!]), at large scale and long time, behaves as a reaction-diffusion 
equation (i.e. with u = 0) where Dq is replaced by D c g. Therefore for FKPP nonlinearity, 



using Eq. (0) one has: Vf ~ 2y -^2. For a detailed derivation of this statement, the reader 
is referred to the Appendix [A|. 

III. A DISCRETE-TIME APPROACH 

Let us now briefly discuss the general idea of our numerical approach to the study of front 
dynamics in terms of discrete-time maps. 

The physical meaning of Eq. fl5|) is made clear by the limit D Q = 0. In that case, introducing 
the Lagrangian time derivative 

D d „ 
— = — + u • V . 

Dt dt 



Eq. flip reduces to 



Denoting by F* the formal evolution operator of Eq. (|j) without noise (Dq = 0), i.e. x(£) = 
F*x(0), and by G l the evolution operator of ^8 = ^f(0), i.e. 6{t) = G t 9(0), one can write 
the solution of Eq. @ in the form: 

0(x,t) = G*0(F-*x,O) . (10) 

Equation ( ]TD| ) is nothing but Eq. (^) in the absence of molecular diffusivity, i.e. when only 
one path ends in x at time t. 
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In the following, we will concentrate on laminar velocity fields and we will develop a suitable 
framework to compute some essential properties for these systems. In time-periodic velocity 
fields, u(x, t+At) = u(x, t) where At is the period and for Do = the Lagrangian motion can 
be described by a discrete-time dynamical system. In other words, the position x(t + At) 
is univoquely determined by x(t); in addition, because of the time-periodicity, the map 
x(t) — > x(t + At) does not depend on t. We remind that a periodic time dependence is 



sufficient to induce Lagrangian chaos [p4 



With these considerations in mind, we can write a Lagrangian map for the position: 

x(t + At) = Fa*(x(*)) . (11) 



If u is incompressible, the map ( |TT]) is volume preserving i.e. its Jacobian matrix has unit 
determinant. In the following we will limit our analysis to the bidimensional case. In that 
situation the map ( |TTD is symplectic, and the dynamics is described by a discrete-time version 
of an Hamiltonian system. Of course it is not simple at all to find explicitly Fa*(x) for a 
generic velocity field. On the contrary, it is not difficult to build Fa* in such a way that the 
qualitative features of a given flow are well modeled, as we will show below. 
Another situation in which one obtains exactly a discrete-time map (fll|) for the Lagrangian 
motion, is the case of velocity field which is always zero apart from 5-impulses at times 
t = 0, ±At, ±2 At, ±3 At, . . . 



u(x,t)= ^ v(x)6(t - nAt) . (12) 

n=— oo 

The effects of a non zero diffusivity are taken into account by adding a noise term 

x(t + At) = F At (x(t)) + ^2D At w(t) , (13) 

where w(t) are standard independent gaussian variables. 
If the production term also is zero apart from (^-impulses 

oo 

f(0)= 9(0)S(t-nAt) , (14) 

n=— oo 
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one can introduce a reaction map 

6(t + At)=G At (6(t)) . (15) 

Now we are ready to write the dependence of the field 9 at time t + At on the field at time 
t in terms of the advection and reaction maps, F At and G At 

0(x, t + At) = (G A ^(F A J(x - y/2D At w(t)), t))^ (16) 

Equation ( |16D is exactly equivalent to Eq. fl5[) for maps (for velocity field and reaction given 
by periodic 5- impulses). 

The concentration field just after the kick, #(x, t + 0), can be written as 

0(x,* + O) = Gf At (0(F^(x),t)) . (17) 

The concentration field #(x, £ + At — 0) is obtained from #(x, t + 0) solving the bare diffusion 
equation d t 9 = D V 2 9 with the initial condition given by Eq. (|T7|): 

0(x,t + At-O) = y e-^0(x - y/2D At w, t + 0) dw , (18) 

which is nothing but Eq. flT6|), and d is the dimension of the space, x e H d . 
From an algorithmic point of view the whole process between t and t + At can be thus 
divided into three steps, a diffusive, an advective and a reactive one. The first two steps 
determine the origin of the Lagrangian trajectory ending in x and accordingly have to evolve 
backwards in time with a given noise realization w. In the third step, the reaction at the 
point x for the advected/diffused passive scalar 9 is computed: 

1) backward diffusion: x — ► x — w 

2) backward advection via the Lagrangian map: x — y/2D At w — >• F _1 (x — y/2D At w) 

3) forward reaction: 9{t + At) = G At (9(t)). 

Let us remark that Eq. ( |TED is exact if both the velocity field and the reaction are 5-pulsed 
processes. However one can also use the formula ( |l6l) as a practical method for the numerical 
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integration of Eq.fjj]) if one assumes small enough At, so that the Lagrangian and reaction 
maps are given at the lowest order by 

F At (x) ~x + u(x)At , G At ~e + —f(9). 

A. The choice of the reaction map 

We now introduce a reaction map G A t(9) corresponding to FKPP nonlinearity. The map is 
characterized by an unstable fixed point in 9 = 0, a stable one in 9 = 1, and a convex shape: 

G A t{9) = (1 + aAt)9 + 0(At 9 2 ) for 9 ~ 

G At {9) = 1 - (3At{l -9) + 0(At(l - 9) 2 ) for 9 ~ 1 

Similarly, for the Arrhenius case we define 

G A t{9) = 9 + 0{At e~ 9c/d ) for ~ 

G At {9) = 1 - /3At(l - 0) + 0(At(l - 0) 2 ) for 9 ~ 1 , 

(see Figure 0). We expect from known results [||], for the time-continuous PDE ([!]) that at 
a qualitative level, the details in the shape of G At {9) are not very relevant, within a given 
class of nonlinearities (e.g. FKPP). This expectation is confirmed by numerical simulations. 
Naturally, if one is interested in the details of some specific combustion processes, one has 
to work with a precise shape of G At . 

B. The choice of the lagrangian map 

If we limit our study to the 2D case, the incompressibility of the velocity field implies 
symplecticity of the map (|TT|). A rather general class of symplectic maps is the following: 



x{t + At) = x{t) + PAt {y{t)) 

y(t + At) = y(t) + q At (x(t + At)) . 
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It is easy to verify that ( |19"D is symplectic for any choice of PAt(-) and qAt( 4 )- If PAt(y) = 
k At sin?/ and qAt( x ) — k At sinx one has the so called Harper model (often studied in the 
quantum chaos context) corresponding to a chaotic transport in spatially periodic cellular 
structures. The case qAt{x) = gives a non-chaotic shear flow in the x direction. The 
celebrated standard map ftBfl , which is a paradigmatic model for chaotic behaviour in low 
dimensional Hamiltonian system, is obtained with pAt{y) = k At siny and qAt{x) = x. 
In the following we will study two limit cases which give nonchaotic transport, i.e. the role 
of the molecular diffusivity is very important: 

a) open flow field (shear flow) where all the streamlines are open 

b) convective rolls where all the streamlines are closed. 



For the shear flow, we set qAt(%) = in (19), as mentioned above. For the cellular flows we 
have built the Lagrangian map in the following way: consider a 2d incompressible steady 
velocity field u(x,y) = (—d y tfj(x,y),d x ip(x,y)) generated by the streamfunction 

^(*,y) = f/Lsin(^)sin(^) , (20) 

with L-periodic conditions in y and infinite extent along the x-axis. The Lagrangian map 
Fa*(x) is given by the exact integration of the equation ^x = u(x) on an interval At. The 
shape of Fa*(x), i.e. the expression of x(t + At) as a function of x(t) is found explicitly in 
terms of elliptic functions. 

In addition to the cases a) and b) we will study the relevance of a "transversal" perturbation 
to the shear flow. 

IV. NUMERICAL RESULTS 

Since we are interested in the front propagation in one direction, say the x-direction, we 
applied in our simulations periodic boundary conditions in y-direction: 

6(x,y,t) = 6(x,y + L y ,t) (21) 
13 



whereas in x-direction we have 



lim 6{x,y,t) = , (22) 

x— >oo 

lim 6{x,y,t) = 1 . (23) 



a;— >— oo 



In this way the front propagates from left to right. 
The instantaneous front speed V(t) is defined by 

V{t) = J J dx dy[9(x, y, t + At) - 6{x, y, t)\ . (24) 

Since V(t) in general shows some oscillations in time, one can define the mean front speed 
Vf by the time average over a sufficiently long time, after the transient. 
The numerical implementation of the formula (|T6|) is described in detail in Appendix C. 
We have first checked the numerical code to confirm the known results for the front speed 
Vf in the case of FKPP nonlinearity: 

i) only molecular diffusion is present, i.e. F(x) = x. In this case, equivalent to u = 0, 
one has the discrete time version of the FKPP formula (^|) 



Vf = XtV DoAtln{G ' Atm ■ (25) 

See Appendix B for its derivation, 
ii) the reaction is very slow, i.e., r r = At/ ln(G" A< (0)) >> t a (where t a is the advec- 



tion time). In this case homogenization techniques can be applied (see e.g. [17|, and 
Appendix [AT) and one finds 



V f ~ ^D cS At In(G' At (0)) = 2J^ . (26) 



The effective Diffusion coefficient D e ^ in the x-direction is defined by 

D cs = lim (27) 

t^oo Zt 
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and in general must be computed numerically iterating the map ([13]). 

In order to show the validity of the homogenization limit fl26| ) we consider a system where 

the lagrangian motion is given by the standard map 



x{t + l) = x(t) + K sin(y(t)) 
y(t+l) = y(t) + x(t+l) mod 2tt 




(28) 



and the reaction map G(9) is in the FKPP class 

G{6) = 9 + c0(l -9) if c < 1 . (29) 

In the case of c > 1, in order to avoid unbounded reaction map, we use 

- 9) for 9 < 9* 
G{9) = { " (30) 

9) for 9 > 9* 



where 9* = (c + 2 - Vc 2 + 4)/2c and a = (c + 2 - Vc 2 + 4)/(c - 2 + V4 + c 2 ). 

By varying c one can change the ratio between the advection time t a ~ 0(1) and the reaction 

time r r = ln J/( ) = i n (i +c ) ) i- e - the Damkohler number. Figure |2| shows Vf vs 1/c, for K = 1 

and K = 3. At large r r , i.e. large c -1 the homogenization limit is recovered. 

We now move to the description of the numerical results in three case studies. 



A. Shear Flow 

The shear flow is the simplest case to study and it will be presented here shortly. For our 
simulations we use the system (|1J) with q = 0: 

x{t + 1) = x{t) + U sin(27n/(t) / L y ) 
y(t + l)=y(t) . 

and a reaction map G{9) of the shape ( p9|) (p0|), where for the sake of simplicity we use 
At = 1. While keeping fixed the diffusivity we investigate the different regimes in the (U, r r ) 
space. Among all the possible combinations of diffusion, advection and reaction time scales, 



15 



we assume in this paper always that diffusion is the slowest one. 



The front velocity in the homogenization regime, i.e. very slow reaction, is Vf — 2y 
where D e g can be easily computed for the shear flow |30[]: D^ — Dq ~ U 2 /D . Therefore for 
slow reaction one has a linear behaviour Vf ~ U. The bounds discussed in the introduction 
suggest that in general Vf = all + b where a and b may depend on r r and D , as confirmed 
by recent numerical results |3l] . 

In Fig. H] we show two snapshots of the concentration field for slow and fast reaction. In 
Fig. f| the front velocity is displayed in dependence on the advection velocity U for different 
reaction times, homogeneization holds for slow reaction rates; decreasing r r the front speed 
increases until for high reaction rates the geometrical optics regime is reached. 
Our results obtained with the discrete time map approach are in perfect agreement with 
the direct numerical simulations presented in Ref. |31| , where the dynamical equations are 
solved in the Eulerian framework. 



B. Cellular Flow 

The numerical simulations have been performed using the velocity field u(x, y) = 
(—dyip(x,y),d x ip(x,y)) generated by the streamfunction defined in Eq. fl20l) and G{9) given 
by (|29|) and (|30|) with c = At/r r and At <C r r . In Fig. |] we show a snapshot of the 
concentration field for two values of r r . 

The key to the understanding of the different regimes of front propagation stands in the 
description of front dynamics in terms of effective macroscopic equations, which we introduce 
hereafter, following Ref. ||36|| . The dynamics of 9 is characterized by the length-scale of the 
cell-size, L. We can therefore perform a space discretization which reduces each cell, Ci, to 
a point, i, mapping the domain - a two-dimensional infinite strip - to a one- dimensional 
lattice, and the field 9 to a function defined on the lattice 0, = L~ 2 j c 9 dxdy. Integrating 
Eq. ([1]) over the cell Ci we obtain Oj = Jj + i — J» + \% where Jj = L~ 2 f, ft D d x 9 dy is the 
flux of matter through the left boundary of the 2-th cell, and Xi = L~ 2 f c r~ l f{9) dxdy is 
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the rate of change of 0j due to reaction taking place within the cell. We will show that 
it is possible to model the dynamics with a space-discretized macroscopic reaction-diffusion 
equation: 

e, = D eS (~e i+1 - e, + -e^) + ^(e.) . (32) 

The effect of velocity is to renormalize the values of diffusivity, -Do ~~ ► D c g(D ,U, L), and 
reaction time-scale, r r — > r e ff(r r , U, L) and therefore the advective term does not appear 
any longer in the effective dynamics given by Eq. (|32"|). The assumption that k and r are 
independently renormalized by advection is consistent in the regime v/U = (Da/ Pe) 1 / 2 <C 1 
The renormalized diffusivity D e g accounts for the process of diffusion from cell to cell result- 
ing from the nontrivial interaction of advection and molecular diffusion. The renormalized 
reaction-time r e g amounts to the time that it takes to a single cell to be completely burned, 
and depends on the interaction of advection and production. In that context, the limiting 



speed of the front in the moving medium will be Vf ~ 2^J D c ^/t c ^. The goal is now to derive 
the expressions for the renormalized parameters from physical considerations. 
Renormalization of diffusivity. To obtain the value of D e g it is sufficient to neglect the 
reaction term in equation (ffl), i.e. consider a passive scalar in a cellular flow. The solution 



is known, fl2|-|l4| 



^ ~ Pe 1 ' 2 Pe » 1 . (33) 
Do 

For large Pe (D small) the cell-to-cell diffusion mechanism can be qualitatively understood 
in the following way: the probability for a particle of the scalar to jump across the boundary 
of the cell, in a circulation time L/U, by virtue of molecular diffusion, can be estimated as 



the ratio of the diffusive motion across the streamlines, 0(WDqL/U), to advective motion 
along streamlines, O(L), leading to p ~ (D /(t/D)) 1//2 , hence the effective diffusivity _D eff ~ 
pUL~ D Pe 1/2 . 

Renormalization of reaction time. At small Da, where reaction is significantly slower than 
advection, the cell is first invaded by a mixture of reactants and products (with a low content 
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of products, Gj 1) on the fast advective time-scale, and complete reaction (O, = 1) is 
then achieved on the slow time-scale r c g ~ r r (Figure |6|). The area where the reaction takes 
place extends over several cells, i.e. the front is "distributed" . 

At large Da, the ratio of time-scales reverses, and in a (now short) time r r two well-separated 
phases emerge inside the cell. The interface has a depth A ~ \JD§t t = LPe~ l / 2 Da~ 1 / 2 i.e. 
it is thin compared to the cell size. Here the process is characterized by an inward spiral 
motion of the outer, stable phase (see Figure 0), at a speed proportional to U as it usually 
happens for a front in a shear flow at large Da. Indeed it is easy to show that, inside a 
cell, the problem can be mapped to a front propagation in a shear-flow in "action-angle" 
variables Therefore the 6 = 1 phase fills the whole cell on the advective time-scale, 

giving t cS ~ L/U. 

In summary, we have the following behavior for the renormalized reaction time 



Tcff 



1 Da < 1 

(34) 

Da Da > 1 . 



Now, we have all the information to derive the effective speed of front propagation for a 



cellular flow. Recalling that Vf ~ \/ D c s/r e s, we have for the front velocity the final result 

V f Pe 1/A Da « l,Pe> 1 

(35) 

Pe 1/4 Da' 1/2 Da > 1, Pe > 1 



where we restricted ourselves to the most interesting case Pe 3> 1. At small Da the front 
propagates with an effective velocity which scales as the upper bound derived above, that is 
as Pe 1 / 4 . At large Da front speed is less enhanced than at small Da: according to Eq. (j35|). 
we have Vf/ ^4D c g/T r ~ Da^ 1 / 2 for Da ^> 1. In terms of the typical velocity of the cellular 
flow, we have Vf oc U 1 ^ for "fast" advection (U 3> L/r r , or equivalently Da <^ 1) whereas 
Vf oc U 3 / 4 for "slow" advection (U <C L/r r , or Da 3> 1). The numerical results are shown in 
Figure || The case of "fast" advection corresponds to the one with slow reaction, for which 
the homogenization limit holds. 
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In the geometrical optics limit Da 3> l,Pe ^> 1, the effective speed of the front is pro- 
portional to the area of the interface that separates the two phases. In two dimensions, 
the interface is characterized by its length, £, and its depth, A. We have the relationship 
Vf ~ X£/(Lr r ) which entails the result that the ratio of the length of the interface in a 
moving medium, £, to the length in a medium at rest, L, is £/L ~ Vf/Vo and is therefore 
larger than unity. The structure responsible for this elongation of the front edge is the spiral 
wave shown in Figure [7[ 

Finally, it is interesting to look at the shape of the effective reaction term r~gF(Q) appearing 
in the renormalized equation (0). As shown in Figure for small Da the effective produc- 
tion term is indistinguishable from the "bare" one. Increasing Da, the reaction rate tends 
to reduce progressively, inducing the slow-down of the front speed. The effective potential 
shows a small region where the production term is essentially the microscopic one, followed 
by an intermediate regime characterized by a linear dependence on the cell-averaged con- 
centration, with a slope directly proportional to Da~ x . That is in agreement with a typical 
effective reaction time r e fj ~ r r Da (cf. Eq. (P^)). 

We conclude the discussion on cellular flows by noting that the scaling behaviors Vf vs U 
are in agreement with the rigorous bounds Vf > C\U 1 ^ + C 2 and Vf > C3U 1 / 2 + C 4 for slow 
and fast reaction |18|| . 



In the previous subsections we discussed pure cellular and pure shear flows. Now we in- 
vestigate the transition between these two limiting cases. To this aim, we will use for the 
lagrangian motion the generalized Harper map: 



The case Ut = corresponds to shear flow, whereas Ut = U gives a chaotic cellular flow. 



C. Percolating Flow 




y(t + 1) = y(t) + U T smx(t + 1) mod 2vr . 



(36) 
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In order to give an idea of the Lagrangian behaviour of the flow generated by ( |36| ) we show in 
Fig. some trajectories at different value of Ut- For Ut 7^ the map ( |36| ) exhibits chaotic 
behaviour in some regions. At small value of Ut one has basically a ballistic transport in the 
x direction apart from small recirculation regions. For Ut ~ U a typical "cat's eye" pattern 
appears with percolating channels among the re-circulation regions. A chaotic cellular flow, 
rather similar to the case of convective cells discussed in section IV B (apart from a rotation 



of 7r/4) is obtained for Ut = U, see Fig. [TO). The behaviour for Ut 3> U can be understood 
by a simple statistical argument valid for large K |38|| : at very large values of Ut, y{t) 
changes very rapidly, therefore the oscillatory term sin(?/(£)) can be considered as a zero 
mean random process and the variable x(t) is well approximated by a diffusive process with 
D eS ~ U 2 /A. 

In Fig. [11] we show snapshots of the concentration field for different reaction times. The 
natural question is how the transition from pure shear to percolation and from cellular flow 
to percolation, respectively, changes the front speed. On an intuitive basis we expect the 
front propagation to be slower in the cellular case than in the shear case. In Fig. |T2 



we 



plot the front speed in dependence of the sidewind Ut for different reaction time scales. 
The diffusion coefficient of the system (RBI) has been calculated in separate runs for the 
comparison with the homogeneization expectations. 

Let us discuss the figure going from left to right, augmenting the sidewind. For zero sidewind, 
we recover the pure shear result, at the value Ut = U, we recover the pure cellular flow result 
and if we go even beyond, for a very large sidewind, the reaction becomes relatively small 
and thus we enter the homogeneization regime. This transition pure shear-pure cells is 
smooth as confirmed by the figure. The three horizontal lines show the asymptotic values of 



2A/-D eff /r r (the value of D eS does not change significantly as a function of Ut if Ut is large 
enough) . 
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D. Final remarks 



In the previous subsections we have shown the results for the behaviour of Vf as function 
of U and r r for different laminar flows discussing in detail the FKPP nonlinearity. It is 
natural to ask about the effects of the shape of f(9) on the front speed Vf. In particular it 
is interesting to know whether the choice of an Arrhenius nonlinearity changes significantly 
the scenario presented in the previous sections. 

It is known that for ignition nonlinearity, i.e. f{9) = for 9 < 9 C , and expectedly also in the 
Arrhenius case, the flow can suppress front propagation. This effect, called flame quenching, 
is absent for FKKP production terms. This observation may lead to the assumption that 
the front evolution could depend on the shape of f(9) in a dramatic way. However, flame 
quenching takes place only if initial conditions of the field 9 are localized, i.e. 9 is different 
from zero only in a bounded region ||. For the initial and boundary conditions that we use 
here (Eqs. fl2~2|)), the front propagates always fl7,18|, also in the Arrhenius case. For that 
reason, we do not expect major differences in the scaling properties of propagation speeds. 
Indeed, in the particular geometry we use, i.e. an open flow with an infinite reservoir of 
burned material, the values of Vf in the case of Arrhenius nonlinearity are very similar to 
the ones obtained using the FKPP nonlinearity. In Fig. [13] we show Vf as a function of U 
for the cellular flow introduced in Section IV B in the case of Arrhenius nonlinearity. The 
scaling laws Vf ~ U 1 ^ and Vf ~ U 3 ^ for slow and fast reaction hold also in this case. Also 
for shear and percolating flow we do not observe qualitative changes when varying the shape 
of f(9). 

Although the qualitative behaviour of Vf as a function of the system parameters does not 
change for different reaction terms, there are differences in the front speed relaxation to its 
asymptotic value. In the case of FKPP nonlinearity (i.e. "pulled" fronts), without advection 
(u = 0), it is known || that the front velocity relaxes algebraically slow to its asymptotic 
value. Therefore one can expect some numerical difficulties to find out the value of Vf, in 
particular for the slow reaction case in which the front may interest a very large spatial 
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region. With this "caveat" in mind, in our simulations we varied the system size to carefully 
check the convergence of Vf. Using Arrhenius nonlinearity (i.e. "pushed" front) in the case 
without advection is known that the convergence to the asymptotic value is exponentially 
fast. Also in presence of a velocity field we observe that the convergence is much faster than 
in the FKPP case. 

V. SUMMARY AND CONCLUSIONS 

Enhancement of front propagation by an underlying flow is a generic phenomenon for 
advection-reaction-diffusion systems. A relevant question is how the front speed Vf de- 
pends on the detailed properties of the advecting velocity field, in particular on the typical 
velocity U. For an arbitrary flow, it is extremely difficult to derive this dependence analyti- 
cally. Here, we have shown that for all incompressible flows there exists an upper bound to 
the front speed that links it to a single global property of the flow, its effective diffusion coef- 
ficient. The analytic derivation of the effective diffusivity for a given velocity field is itself a 
daunting task, but several generic properties are known and some exact results are available 
for simple flows. In the special case of fast advection, as compared to the reaction timescale 
r r , the upper bound is sharp, and therefore it is possible to obtain the dependence of Vf on 
U . When molecular diffusivity is small we have Vf ~ U for flows with open streamlines such 
as the shear flow, and Vf ~ U 1 ^ for cellular flows. For slow advection, the bound ceases 
to be effective, and one has to resort to numerical simulations in order to determine the 
front speed. We find that for open-streamline flows there is still a linear dependence Vf ~ U 
whereas cellular flows display a Vf ~ U 3 ^ dependence. Within the class of initial/boundary 
conditions for which no flame quenching effect ever takes place, those scaling laws appear 
to be universal with respect to the details of the reaction mechanism. 

Which lessons can we draw from the present results for the open, challenging problem of 
front propagation in turbulent flows? The main point that we want to emphasize is the 
central role played by the effective diffusion process in determining the front speed. That 
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is a reflection of the deep-rooted link between front propagation and transport properties. 
The knowledge of turbulent transport has experienced a significant progress in the past few 



years (see e.g. Refs and PD| ). We believe that those results will reveal helpful to shed 
light on the issue of front propagation in turbulent flow (for a similar point of view, see Ref. 

El)- 
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APPENDIX A: HOMOGENIZATION REGIME 



In this appendix we present briefly the application of homogenization techniques |42|] to the 
reaction-advection-diffusion equation ([I]). We are interested in the large-time, large-scale 
asymptotics for slow reaction timescales. Introducing the small parameter e, we consider 
reaction times r r = f r e~ 2 , and look at the solutions of Eq. ([[]) for times 0(e~ 2 ) and scales 
0(e _1 ). The separation in time-scales allows a multi-scale treatment. Slow variables X = ex 
and T = e 2 t are introduced along with the fast variables x and t. Slow and fast variables 
are considered as being independent. As a consequence space-derivatives act as d x + edx 
and the time-derivative as d t + e 2 dx- The concentration field is expanded in a power series 
in e as 6*(x, t, X, T) = 6*^(x, t, X, T) + e^^x, t, X, T) + . . ., and this expression is plugged 
in Eq. (|l|). At zero-th order in e the equation reads 

dtO® + u ■ V^(°) = D V^W . 
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Due to the dissipative nature of the latter equation, the solution at zero-th order will decay 
to its average value on fast timescales: 

0W(x,t,X,T) =0<°>(X,T) . 

At order e, we obtain the linear equation 

d t e^ + u • v^ (1) - D v 2 x e^ = -u • Vx^ (0) 

which allows the solution 

(x, t, X, T) = 0« (X, T) + w(x, t) • V X (O) (X, T) , 

provided that the auxiliary field w obeys the equation 

d t w + u ■ V x w = _D V 2 .w — u . 

Remark that the production term has not shown up yet. It is at order e 2 that it enters the 
scene 

dt9 m + u • v,0 (2) - DoVle® = -d T e^ - u • v x e {1) + D v 2 x e^+ 

2D V x -V x e^ + -f(9^). 
The solvability condition for the equation at second order requires that 

d T e^ = d v 2 x 9® - (u ■ Vx^ (1) ) + -f(e {0) ) 

T r 

where the brackets denote the average over the fast variables. Plugging the expression for 
into the solvability condition yields the effective, large-scale equation 



where D e s is in general a tensor with components _D c g = DqS^ — h(u l w^ +u^w l ) . Considering 
the propagation in the x-direction only, we recover an effective equation which is equivalent 
to Eq. (|I]) with an effective diffusivity. Therefore the front propagates at a maximal speed 



(computed in the original, fast variables) given by Vf < 2w^2 {sv^o^p-^, for FKPP 



nonlinearity one has Vf = 2 y ^2 ( where _D e ff = -Djj 1 . 
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APPENDIX B: FRONT SPEED FOR DISCRETE TIME MAPS 



The front speed (p5|) for discrete maps can be obtained by simple considerations just following 
the standard way used for the derivation of Vf in the continuous time limit: 

^e = D o v a e + -f{0) . (bi) 

at r r 

Let us consider a front propagating from left to right. For the sake of simplicity we discuss 
a FKPP nonlinearity for the one dimensional case. For x — ► oo, 9(x,t) has an exponential 
shape, 

6{x,t) = e at ~ bx , (B2) 



up to exponentially subleading terms. Inserting ( |B2| ) in ( |Bip and linearizing around 9 = 
one has: 



a = D b 2 + -f'(0) . (B3) 



A saddle point argument gives a selection criterion which allows for the determination of 
the front speed 0: 

^ = min ^ = 2^°/'(0). (B4) 
Consider now the discrete time reaction case, i.e.: 

-0 = D V 2 9 + 9{d)5{t - n) , (B5) 

n=— 00 

where for sake of simplicity we adopt At = 1 (see Eq. (|14])). 
Indicating with G(9) the reaction map (see Eq. (^)) one has: 

9(x,t + 0) = G(9(x,t-0)) . 

Integrating the diffusion equation d t 9 = D V 2 9 between t + and t + 1 — one has: 

9{x, t + 1 - 0) = J e~^G (9{x - ^/2I\ w, t + 0)) dw , (B6) 
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Assuming the shape flB2"|) and linearizing around 9 = 0, i.e.: G{9) ~ G'{0)9, a simple 
gaussian integration gives: 

e a(t+l)-bx ^ e lnG' (0)+D b 2 -bx+at 

The above result implies 

a = lnG'{0) + D b 2 

this is nothing but Eq. (|F33| ) now with In G'(0) instead of ^/'(O). The same selection criterion 
gives 

V f = 2 v / J D lnG / (0) • (B7) 

APPENDIX C: NUMERICAL METHOD 

Since we are interested in propagation along the x-axis, we consider a slab with sides L x ^> 
L y . The boundary conditions are periodic in the ^-direction 9(x,y,t) = 9(x,y + L y ,t). To 
fulfill the conditions (|22|) lim^oo 9(x, y, t) = and lim :I .^_ 00 #(x, y,t) = 1 numerically, we set 
8(0, y, t) = 1 and ^L^, y, t) = 0, which is a good approximation as long as the fronts leading 
edge has not reached L x . We introduce a lattice of step size Ax and Ay (for sake of simplicity 
we assume Ax = Ay) so that the field 8(x, y, t) is defined on the points x n m = (nAx, mAy). 
The numerical code computes 9 n ^ m (t + At) = 9{nAx, mAy, t + At) in terms of 9 n)m {t) using 
Eq. (^). For each grid point x n m , one introduces N independent standard gaussian variables 
W a , a = 1, . . . , N, iV 3> 1, and computes x" m = x n/m — y/2D At W° and finally from this 
r" m = F _1 (x" m ). For 6 , n>m (t + Ai) one needs the values of 6 at time t in the positions r" m . 
In general the n are not on the grid points (nAx, mAy), nevertheless we can compute the 
value 9(r* m ,t) using linear interpolation from 9 nim (t). Therefore we have 

1 - 

9 n>m (t + At) = -J2 G l d Km,t)}. (CI) 

a=l 
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Typically one has a good convergence for N = 50. To simulate the diffusion process we 
have to impose a relation between D , Ax and At to insure that the diffusion transports a 



particle over distances larger than the grid-size, v/2 ?° At > 1 (see Fig III) 
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FIG. 1. Two examples of reaction maps. On the left FKPP type, on the right Arrhenius type. 
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FIG. 2. Front speed for the standard map (28) with reaction map given by Eq. (|29| ) (pOf) , as 
function of c, Dq = 0.04. The upper curve is for K = 3.0, the lower for K = 1.0. The dotted lines 



are the homogenization curves 2y'^) e gdn(T0 + c). The diffusion coefficient D e s depends on K and 
has been computed numerically. 
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FIG. 3. Snapshots of the field 9(x,y) for the shear map (31) and the reaction map (29,30). 
U = 0.5, Do = 0.01, c = 0.2 and c = 2.0 for the upper and lower image respectively. The system 
size is L y = 2ir and L x = 207r. 
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FIG. 5. Snapshots of the field 9(x, y) for cellular flow given by the streamfunction ( ^0|) and the 
reaction map (|29| , |30|) , U = 2.0, Dq = 0.01, r r = 5.0 and r r = 0.5 for the upper and lower image 




FIG. 6. Cellular flow: six snapshots of the field 9 within the same cell, at six successive times 
with a delay r/6 (from left to right, top to bottom), as a result of the numerical integration of 
equation ([[]). Here Da ~ 0.4, Pe ~ 315. Black stands for 9 = 1, white for 9 = 0. 
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FIG. 7. Cellular flow: six snapshots of the field 9 within the same cell, at six successive times 
with a delay (L/U)/6 (left to right, top to bottom). Here Da = 4, Pe = 315. A spiral wave invades 
the interior of the cell, with a speed comparable to U. 
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FIG. 8. Cellular flow: the front speed Vf as a function of U, the typical flow velocity, with 
Do = 0.04. The lower curve shows data at r r = 20.0 (fast advection). The upper curve shows data 
at T r = 0.2 (slow advection). For comparison, the scalings C/ 1 / 4 and U 3 ^ are shown as dotted and 
dashed lines respectively. The horizontal line indicate Vq (the front speed without advection, i.e. 
U = 0) for T r = 0.2. 
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FIG. 9. Cellular flow: the renormalized reaction term T r /r e ^F{Q) vs for three different 
parameters: Da ~ 4 (□), Da ~ 2 (o) and Da ~ 0.4 (x). The continuous line is f{6). The dotted 
and dash-dotted lines have the slopes (0.2 and 0.4) proportional to Da^ 1 in the region of slow 
advection. 
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FIG. 10. Lagrangian dynamics of test particles evolving according to the Harper's map (36) 
for different values of Ut, with U = 1.5. From top-left moving clockwise: Ut = 0.2, Ut = 0.8, 
U T = 1.5 and C/ T = 3.0. 
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FIG. 11. Snapshots of the field 8(x, y) for percolating flow given by Eq. (p6|). The reaction map 
given by (H©. U = 1.5, U T = 0.8, D = 0.01, c = 0.2 and c = 2.0 for the upper and lower 
image respectively. The system size is L y = 2tt and L x = 20tt. 
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FIG. 12. Vf vs Ut, the Lagrangian map is given by (^), the reaction map G{9) is given by 
(|2^ , |30|) . The three curves correspond to c = 0.1, 1.0, 10.0 (from bottom to top), the "horizontal" 
velocity has been fixed to U = 1.5 for all curves. The asymptotic value of 2^/ D c ff ln(l + c) (for 
the horizontal direction) is shown by the three horizontal lines, the corresponding value of D e $ has 
been calculated numerically for large Ut- 
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FIG. 13. Cellular flow: the front speed Vf as a function of U, the typical flow velocity in the 
case of Arrhenius production type: f(9) = (1 — 9) exp(— 9 c /0). The lower curve shows data at 
r r = 2.0 and 9 C = 0.5 (fast advection). The intermediate curve at r r = 2.0 and 9 C = 0.2 shows 
the crossover from fast (right side) to slow advection (left side). The upper curve shows data at 
r r = 0.2 and 9 C = 0.2 (slow advection). For comparison, the scalings U 1 / 4 and U 3 / 4 are shown as 
dotted and dashed lines respectively. 
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FIG. 14. Pictorial scheme of the numerical algorithm, as discussed in Appendix C. Here, 
if- = y/2Do At W a where W a is a standard gaussian variable. 
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